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The impacts of global climate change on different aspects of 
humanity's diverse life-support systems are complex and often 
difficult to predict. To facilitate policy decisions on mitigation 
and adaptation strategies, it is necessary to understand, quantify, 
and synthesize these climate-change impacts, taking into account 
their uncertainties. Crucial to these decisions is an understanding 
of how impacts in different sectors overlap, as overlapping 
impacts increase exposure, lead to interactions of impacts, and 
are likely to raise adaptation pressure. As a first step we develop 
herein a framework to study coinciding impacts and identify re- 
gional exposure hotspots. This framework can then be used as 
a starting point for regional case studies on vulnerability and mul- 
tifaceted adaptation strategies. We consider impacts related to 
water, agriculture, ecosystems, and malaria at different levels of 
global warming. Multisectoral overlap starts to be seen robustly at 
a mean global warming of 3 °C above the 1980-2010 mean, with 
11% of the world population subject to severe impacts in at least 
two of the four impact sectors at 4 °C. Despite these general con- 
clusions, we find that uncertainty arising from the impact models 
is considerable, and larger than that from the climate models. In 
a low probability-high impact worst-case assessment, almost the 
whole inhabited world is at risk for multisectoral pressures. Hence, 
there is a pressing need for an increased research effort to develop 
a more comprehensive understanding of impacts, as well as for 
the development of policy measures under existing uncertainty. 


amplify effects, restrict response options, and lead to indirect 
impacts in other regions, thus strongly increasing the challenges 
to adaptation (1). In this article we take an important first step 
toward the analysis of these effects through a consistent assess- 
ment of the geographical coincidence of impacts as multisectoral 
exposure hotspots. The Intersectoral Impact Model Intercom- 
parison Project (ISI-MIP, www.isi-mip.org) offers a unique op- 
portunity for this analysis by providing multimodel ensembles of 
climate-change impacts across different sectors in a consistent 
scenario framework. 

Through the investigation of biophysical impacts of climate 
change, which form the linkage between climate and society (2, 
3), this study moves beyond previous hotspot analyses that have 
mostly used purely climatic indicators (4-7). In addition, the set- 
up enables an assessment of uncertainty because of both multiple 
Global Climate Models (GCMs) and multiple Global Impact 
Models (GIMs) in each sector (8). Finally, impacts are analyzed 
at different levels of global mean temperature (GMT) for a 
comparison at different levels of global warming. This global 
analysis serves two objectives. First, tangible adaptation strate- 
gies require knowledge of local vulnerability, defined by expo- 
sure, sensitivity, and adaptive capacity. The regional exposure 
hotspots can therefore serve as a starting point for prioritized 
case studies and studies of interactions as the basis for the de- 


coinciding pressures | differential climate impacts | ISI-MIP 

O ver the coming decades, climate change is likely to signifi- 
cantly alter human and biological systems, pushing the 
boundaries of variability beyond historic values and leading to 
significant changes to what are considered typical conditions. 
Identifying the locations, timings, and features of these impacts 
for a given level of global warming in advance allows the de- 
velopment of appropriate adaptation strategies, or can motivate 
decisions to mitigate climate change. Although climate-change 
impacts are extensively studied in individual sectors, their over- 
laps and interactions are rarely taken into account. However, 
these impacts are likely to be of great consequence, as they can 


velopment of adaptation strategies that can be expanded to ad- 
ditional regions as needed. Second, the focus on GMT change is 
crucial when studying costs and benefits of mitigation policies, 
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such as the 2 °C target set by the international community to 
reduce risks from climate-change impacts and damages (9, 10). 

The analysis comprises four key impact sectors: water, agri- 
culture, ecosystems, and health. Health is represented by malaria 
(albeit, there is only one example of health impacts of climate 
change), which does have potentially severe economic con- 
sequences (11). As metrics for the four sectors, we select river 
discharge as a measure of water availability, crop yields for four 
major staple crops (wheat, rice, soy, and maize) on currently 
rain-fed and irrigated cropland (12) (Fig. SI), the ecosystem 
change metric T (13), and the length of transmission season 
(LTS) for malaria. Although these four metrics do not cover the 
full range of possible societally relevant climate-change impacts, 
they do include crucial aspects of livelihoods and natural 
resources, especially for developing countries: water availability, 
food security, ecosystem stability, and a key health threat. 

We aim to define levels of change in each sector that can be 
considered severe as basis for multisectoral hotspots. “Severe” is 
taken to mean a shift of average conditions across selected 
thresholds representing significant changes relative to the his- 
torical norm. A multisectoral perspective is thus possible through 
the simultaneous occurrence of above-threshold changes in 
multiple sectors. Although climate change can have both positive 
and negative impacts, for the purpose of vulnerability analysis we 
identify hotspots of changes that put additional stresses on hu- 
man and biological systems. Average conditions are measured as 
the median over 31-y time periods. For the thresholds, we take 
a statistical approach for water availability and crop yields, 
whereas we use a more comprehensive metric for ecosystem 
change, and resort to a relatively simple indicator for malaria 
conditions. The thresholds in the water and agricultural sectors 
are defined as the 10th percentile of the reference period dis- 
tribution (1980-2010) of discharge and crop yields, respectively. 
This threshold means a shift of average conditions into what is 
considered today moderately extreme, happening in only 10% of 
all years. Behavior is robust to the choice of a smaller threshold 
(Fig. S2). This low end of the distribution excludes floods, as the 
focus is on reduced water availability. Clearly, the chance to 
cross the threshold depends on the level of variability in a given 
region and may in fact mean relatively small absolute change; 
however, it reflects the assumption that people in regions already 
subject to highly variable conditions are better prepared to adapt 
to more extreme average conditions (14). 

The T-metric (13) represents the difference between future 
states of ecosystems and present day conditions through an ag- 
gregate measure of changes in stores and fluxes of carbon and 
water, as well as vegetation structures. A large value of T indi- 
cates significant changes in biogeochemical conditions or vege- 
tation structure, which would likely lead to considerable trans- 
formations of the ecosystem. Based on differences between 
present day ecosystems, Heyder et al. (13) define T > 0.3 as the 
threshold for a risk of severe change, [see also Warszawski et al., 
(15)]. Such changes may reduce biodiversity, which is crucial for 
the resilience of many ecosystem services (16). Furthermore, the 
livelihoods of many vulnerable populations, along with cultural 
values and traditions, are closely tied to existing ecosystems (17). 
The threshold for changes in the prevalence of malaria is defined 
as a shift in the LTS, from < 3 mo to >3 mo. This shift corresponds 
approximately to a switch from epidemic to endemic malaria based 
on climatic conditions (based on data from the Mapping Malaria 
Risk in Africa project, www.mara.org.za) (Fig. S3). 

All impacts are simulated with multiple, predominantly pro- 
cess-based GIMs (agriculture and ecosystems, 7 models each; 
water, 11 models; malaria, 4 models). These GIMs are driven by 
three GCMs, simulating the highest representative concentration 
pathway (RCP8.5) (18). Although current emissions are follow- 
ing a similar trajectory, we choose RCP8.5 primarily to cover the 
largest possible temperature range, not as a worst-case scenario 


(19). For each GIM-GCM combination and at each grid point, 
we define a “crossing temperature” that is the GMT change 
(AGMT) at which the sectoral metric crosses the respective 
impact threshold. Sectoral crossing temperatures are then taken 
as the median over all GIM-GCM combinations of a given sector. 
In our strict assessment, only robust results are taken into ac- 
count, defined as an agreement of at least 50% of all GIM-GCM 
combinations of a given sector at which the threshold is crossed. 
Overlapping pressures at a given grid point are assumed to arise 
when multiple sectors have crossed at a given AGMT. Results 
are presented in terms of total area affected by the shift as a 
function of AGMT. Note that GMT changes in this report are 
with respect to the 1980-2010 period, which is ~0.7 °C above 
preindustrial levels (20). 

Results and Discussion 

Sectoral Analysis. The basis for the study of multisectoral overlap 
is the AGMT level at which the thresholds for severe change are 
crossed (if at all) in each of the four sectors (Fig. 1 and Fig. S4). 
Median 31-y water availability is projected to drop below the 
reference distribution’s 10th percentile in the Mediterranean, 
regions of South America, in particular the southern Amazon 
basin, regions in coastal western and central Africa, and parts of 
south-central Asia for a warming of up to 4.5 °C under RCP8.5. 
This distribution includes some regions of large projected rela- 
tive drop in discharge (21), although the relatively strict 10th 
percentile criterion means that it does not capture all of them (e. 
g., southern United States). The regions affected by crop yields 
below the threshold are tropical regions dominated by rain-fed 
agriculture; this is consistent with the expectation that rain-fed 
systems are likely to see larger and more consistent yield losses 
than irrigated areas that can adapt more successfully. No nega- 
tive effects on yields are seen at higher latitudes, as these initially 
benefit from higher temperatures and CO 2 fertilization effects 
and exhibit yield increases (22, 23). For both discharge and 
yields, thresholds start to be crossed at AGMT = 1 °C. 

Significant risk of ecosystem change, as indicated by the 
T-metric, has the largest geographical extent of all sectors, with 
most regions exhibiting crossing temperatures of 3-4 °C. This 
large extent occurs because it encompasses very different eco- 
system responses. There is forest die -back because of less rainfall 
in the Amazon and heat stress in boreal forest regions, but also 
increased greening in Europe and Africa because of warmer, 
wetter conditions, as well as replacement of some vegetation 
species with others better adapted to the new conditions. Forest 
advances northward as a result of higher temperatures and the 
trees’ increased water-use efficiency in response to higher at- 
mospheric CO 2 concentrations. On the Tibetan Plateau, distin- 
guished by the lowest crossing temperature of AGMT = 2 °C, 
increased vegetation growth because of longer growing seasons 
and warmer winters puts the current grass and shrublands at risk. 
Although not all of these changes will be negative per se, they 
would constitute a disruption and possibly a need for adaptation 
of local societies to the prevailing ecosystem conditions. 

Finally, malaria prevalence is expected to increase in higher 
latitudes, higher altitudes, and in regions on the fringes of cur- 
rent malaria regions because of warmer and wetter climatic 
conditions. However, when conditions become drier, prevalence 
can also decrease. As a result of the very different parameter- 
izations used in the four malaria models considered here, 
agreement among models on the changes is poor, leaving very 
few areas as robustly crossing the 3-mo LTS threshold. Never- 
theless, in agreement with previous work, the Ethiopian High- 
lands are one of these regions (24, 25). 

Multisectoral Hotspots. We define hotspots as regions of multi- 
sectoral exposure where two or more of the sectoral metrics have 
crossed their respective thresholds of severe change in average 
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crossing temperatures 

■ AGMT<0.5°C ■ 0.5<AGMT<1.5°C ■ 1.5<AGMT<2.5°C ■ 2.5<AGMT<3.5°C ■ 3.5<AGMT<4.5°C 


Fig. 1. Threshold crossing temperatures with respect to the reference period GMT for the four sectoral metrics: discharge (A), crop yields (6), risk of severe 
ecosystem change (C), and LTS of malaria (D). Areas in white do not cross the respective threshold. The gray color indicates regions which are either masked 
out [discharge, T, crop yields (only regions where the maize, wheat, soy, and rice are currently cultivated are considered)], or where malaria is already endemic 
(D). An agreement of 50% of all GIM-GCM combinations on threshold crossing is required for consideration in the analysis. 


conditions under the strict assessment, which means with high 
likelihood (Fig. 2). According to our results there is no overlap of 
severe change in all four sectors. The most prominent hotspot is 
the southern Amazon basin, with some parts projected to ex- 
perience severe changes in three sectors (yields, ecosystems, and 
discharge) and large areas affected by two pressures. The second 
largest hotspot region is southern Europe, with overlapping 
changes in discharge and ecosystems. These two areas, as well as 
smaller tropical hotspot regions in Central America and Africa, 
were also identified in other studies using different methods, 
supporting our findings (5, 6). In addition, we identify the 
Ethiopian highlands as a hotspot because of the overlap of 
malaria extension, crop yield reduction, and ecosystem change; 
northern regions of south Asia are affected by either reductions 
in discharge and crop yields or crop yield reduction and eco- 
system change. These multisectoral hotspots occur in both regions 
with high population density (i.e., Europe, east Africa, south 



Fig. 2. Multisectoral hotspots of impacts for two (orange) and three (red) 
overlapping sectors in the strict assessment, with 50% of GIM-GCM combi- 
nations agreeing on the threshold crossing in each sector, for a GMT change 
of up to 4.5 °C. Which sectors overlap depends on the location and can be 
discerned from the sectoral patterns in Fig. 1 . An overlap of all four sectors 
does not occur in the strict assessment. Regions in light gray are regions 
where no multisectoral overlap is possible at all because of sectoral restric- 
tions as shown in Fig.1. The dark gray shows the additional regions affected 
by multisectoral pressures under the worst-case assessment, where a mini- 
mum of 10% of all sectoral GIM-GCM combinations have to agree on the 
threshold crossing. 


Asia) and sparsely populated areas (i.e., Amazon). These hot- 
spots cover developed, emerging, and developing economies, 
each with different degrees of adaptive capacity and sensitivity to 
the multisectoral pressures. Note that these factors are not taken 
into account here. A weighting of the relative importance of the 
sectoral pressures depends strongly on local factors, such as so- 
cietal structures and values, economic base, and environmental 
imperatives. Therefore, a more detailed interpretation of the 
hotspots requires in-depth regional case studies, but is beyond 
the scope of this study. 

Regions typically expected as high-exposure regions, like Africa, 
do not emerge strongly as hotspots here, which is partially be- 
cause of the sectors used in the analysis and the individual 
characteristics of the sectoral metrics, both influencing their 
combination. In particular, the global area where three or four 
regions can potentially overlap is limited to where the four staple 
crops are currently cultivated and where malaria is not yet en- 
demic (excluding gray areas in Fig.l). Hence, a different picture 
might arise if, for example, changes in the occurrence of extreme 
events, like droughts and floods, were included as metrics, which 
would likely increase the occurrence of hotspots in Africa and 
south-east Asia (26). 

The Role of Uncertainty. An additional factor limiting the overlap 
of areas with severe change in different sectors is the large un- 
certainty in projections, stemming mainly from the GCMs and 
GIMs. When the results for the three GCMs are separated, 
different multisectoral hotspot patterns emerge, with some re- 
gions only appearing as hotspots with a single GCM (Fig. S5). 
This appearance is because of different sectoral patterns asso- 
ciated with each GCM as a result of variances in projections of 
key climate variables influencing the impact models. Climate 
model uncertainty is therefore an important cause of the limited 
sectoral overlap in our analysis. Uncertainty from impact models, 
however, is much larger (Fig. 3 and Fig. S2). This finding is in 
agreement with previous literature and other analyses in this 
Special Features issue of PNAS (15, 21, 27, 28). Agreement is 
highest among the ecosystem models, whereas differences are 
largest between the global crop models. In addition to un- 
certainty as to whether the thresholds are crossed, there is also 
uncertainty on the crossing temperature, with SDs of around 1 °C 
in most sectors (Fig. S6). The details of the model differences 
are beyond the scope of this report. However, we emphasize the 
importance of accompanying this study with detailed sectoral 
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Fig. 3. Cumulative fraction of global land area (excluding Antarctica; for crop yields the relevant area is the maximum crop area as covered today by the four 
staple crops: maize, wheat, soy, and rice) having crossed the respective sectoral thresholds up to the given AGMT for discharge (A), crop yields (6), risk of 
severe ecosystem changes (C), and LTS (D). Black boxes show the uncertainty across impact models, and red boxes indicate the uncertainty across GCMs. Each 
box indicates the interquartile range, the thick line shows the median, and the whiskers extend over the whole range of the distribution of all GCMs/GIMs at 
that temperature bin. Note the different ranges on the y axis for each panel. 


understanding and analysis, which can be found elsewhere in this 
issue (see also SI Text) (15, 21, 22, 24). 

This high level of uncertainty warrants the strict robustness 
limit of 50% agreement among GIM-GCM combinations used for 
the identification of hotspots. At the same time, this uncertainty 
may mask a remaining risk, given that models appearing at the 
ends of the distribution cannot be disregarded because no per- 
formance-weighting of models was carried out. Therefore, we 
also provide a worst-case assessment of multisectoral hotspots, 
with crossing temperatures determined as the 10th percentile of 
all crossing temperatures in a given grid cell. This process means 
that only 10% of all GIM-GCM combinations have to agree on 
the threshold crossing (chosen to have at least two in a sector, to 


avoid spurious effects of one outlier) and the resulting crossing 
temperatures are lower limits. This worst-case assessment shows 
a large additional extent of multisectoral overlap (Fig. 2, dark 
gray areas) with almost all of the world’s inhabited areas affected. 
The areas with highest exposure in this case have an overlap of all 
four sectors (Fig. 3 and Fig. S7). This worst case is rather extreme, 
but nonetheless it represents the upper end of the risk spectrum in 
light of the large uncertainties. 

Aggregate Effects with GMT. The total global area and population 
that are projected to face average conditions that are considered 
rare today in more than one sector increases with GMT (Fig. 4). 
The likelihood for multisectoral overlap increases with the area 
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Fig. 4. Cumulative fraction of the global area (brightly tinted bars, excluding Antarctica) and population (lightly tinted bars) affected by the thresholds 
being crossed at AGMT in at least two (red), three (orange), and four (blue) overlapping sectors. (Left) The strict case (agreement of at least 50% of GIM-GCM 
combinations on the threshold crossing). (Right) The worst-case (agreement of at least 10% of GIM-GCM combinations on the threshold crossing and the 
sectoral crossing temperature is the 10th percentile of all crossing temperatures) assessment. Overlap of four sectors does not occur in the strict case. 
Population is held constant at the year 2000 levels. 
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affected in the individual sectors, one reason for the onset of 
multisectoral pressures at relatively high levels of GMT change 
only. For the strict assessment, multisectoral severe pressure 
begins at AGMT = 3 °C above the 1980-2010 baseline; at 4 °C 
roughly 6% of the global area (excluding Antarctica) and close to 
11% of the global population are affected. Correspondingly, the 
largest increases in the areas having crossed the thresholds in 
each sector are seen between AGMT = 1 °C and 3 °C (Fig. 3), 
with some indications for saturation after that. Further increases 
in affected areas are possible at higher AGMT levels than 
studied here. For example, a peak and possible later decline in 
crop yields is expected as a result of heat stress overtaking the 
initial benefits of the CO 2 fertilization effect. 

In the worst-case analysis, area and population affected is al- 
ready much larger at lower levels of AGMT, with the largest 
increase between 1 °C and 2 °C above the 1980-2010 baseline 
and an inflection of the trend after that. Almost the entire global 
population is exposed to multisectoral pressure at AGMT = 4 °C. 
In addition, roughly 18% of the global population is projected to 
experience severe pressure in all four sectors. The affected 
regions are in Europe, North America, and south-east Asia (Fig. 
S7), driven by the extension of malaria prevalence to higher 
latitudes. This interpretation may give too much emphasis to this 
pressure, as malaria distribution also depends strongly on so- 
cioeconomic factors but is here only driven by climate suitability 
(29). Nevertheless, the increased overlap of three or even four 
sectors in the worst-case assessment indicates a strong adapta- 
tion pressure, albeit at low probability. 

Implications and Further Research 

This identification of multisectoral hotspots of climate change 
impacts is to our knowledge unique in its use of a consistent 
framework with multiple impact models per sector and using 
AGMT as a metric for climate change. Our global analysis 
provides a starting point for more detailed understanding of the 
extended implications of climate change for exposure and ad- 
aptation actions. Although geographically overlapping impacts 
only start at AGMT = 3 °C above the 1980-2010 baseline (almost 
4 °C above preindustrial GMT levels), large increases in exposed 
areas within the sectors start at around 2.2 °C above preindustrial 
levels. In the worst-case analysis, the largest increase in affected 
area and population occurs between roughly 2 °C and 3 °C above 
preindustrial levels. This finding provides important insight for 
mitigation strategies. 

The identified multisectoral hotspots are geographically di- 
verse, including the southern Amazon basin, southern Europe, 
the Ethiopian highlands, and northern India, and are driven by 
different combinations of coinciding sectors. Implications and 
possible feedbacks between the overlapping sectors can be in- 
vestigated in regional case studies. At the same time, these 
hotspots could affect distant regions through indirect effects, 
such as trade or migration. Appropriate adaptation planning that 
considers coinciding (and also interacting) pressures facilitates 
the development of strategies designed to address such multiple 
challenges, and avoids creating solutions for one pressure that 
possibly seriously exacerbates another (e.g., draining wetlands to 
reduce malaria in an area prone to increases in flooding). 

The set-up for our analysis explicitly includes uncertainty in 
both climate and impact models. This format shows that uncer- 
tainties from both GIMs and GCMs are large, limiting the ro- 
bustness of the conclusions; however, it should not hamper 
action at this point, as some level of uncertainty will always 
be present. In particular the low probability-high impact worst- 
case assessment, which shows a very large extent of multisectoral 
pressures starting at lower temperature changes, provides a 
strong motivation for more detailed impacts research. 

Because it is unique, our analysis is a methodological experi- 
ment, to be refined in the light of experience. Indeed, different 


patterns may emerge if different sectors or absolute magnitudes 
of change are included. A comparison of hotspots generated 
with different methodologies will provide valuable insights into 
impact dynamics. The identification of hotspots of positive 
climate-change impacts would create a more balanced and 
comprehensive picture, but requires different metrics to those 
used here. In addition, although a simple overlap of the different 
sectoral metrics is considered here, the challenge for future 
analyses is also to integrate the interactions between the differ- 
ent sectors and indirect effects over large distances, which may 
alter the spatial pattern of hotspots. Examples are interactions 
between water availability and irrigation or ecosystem services, 
and irrigation and malaria occurrence (28, 30). Furthermore, 
a more comprehensive understanding of human vulnerability 
hotspots requires a thorough analysis, combining highly resolved 
indicators of adaptive capacity and sensitivity (which so far seem 
to be lacking) with biophysical hotspot indicators as measures of 
exposure (2, 3). Nevertheless, our study is an important step to- 
ward a consistent integration of multiple sectors in impacts re- 
search, and identifies the risk of sizable hotspots of multisectoral 
pressures under highly plausible levels of global warming. 

Materials and Methods 

Models and Data. For this analysis, simulations were driven by the three ISI- t 

MIP GCMs that exhibit a AGMT= 4 °C by the end of the 21st century 5 g 

(HadGEM2-ES, MIROC-ESM-CHEM, IPSL-CM5A-LR). To improve statistical f | 

agreement with observations, a bias correction was applied to the climate S lo 

data. This bias constitutes an additional source of uncertainty and reduces => 

the spread of present-day GCM climatologies (31-34). The gridded year 
2000 population data are based on United Nations World Populations 
Prospects data, scaled to match the country totals of the new Shared Socio- 
Economic Pathway population projections for the middle-of-the-road case 
(SSP2; https://secure.iiasa.ac.at/web-apps/ene/SspDb) using the National 
Aeronautics and Space Administration GPWv3 y-2010 (http://sedac.ciesin. 
columbia.edu/data/collection/gpw-v3) gridded population dataset (35, 36). 

Similar results for the percentage of affected global population are found 


spread in results. The global gridded crop model intercomparison was co- 
ordinated by the Agricultural Model Intercomparsion and Improvement 
Project (37). 

Impact Metrics. All metrics have annual temporal resolution, neglecting 
seasonal patterns. To avoid spurious effects, values are set to zero below the 
lower limits 0.01 km 3 yr _1 and 2.5% natural vegetation cover, for discharge 
and ecosystem change, respectively (15, 38). The four crops are combined by 
converting to energy-weighted production per cell using the following 
conversion factors for energy content (MJ kg -1 dry matter): wheat (spring/ 
winter), 1 5.88; rice (paddy), 1 3.47; maize, 1 6.93; soy, 1 5.4 (39, 40). The extent 
of potential agricultural hotspots is limited; for example, millet and sor- 
ghum, which are widely grown in Africa, are not included in the analysis. 
The impact of climate change on malaria occurrence focuses on changes in 
LTS. This simple metric represents an aggregated risk factor because it 
neglects age-dependent immunity acquisition associated with transmission 
intensity. Increases in impacts associated with transitions from malaria-free 
to epidemic conditions are also not considered. 

Hotspots Method. GMT is calculated from the GCM data and change is 
measured with respect to the reference period 1980-2010. The GMT level in 
the reference period is ~0.7 °C above preindustrial, based on estimates for 
1980-1999 of 0.51 °C and the average of the five GCMs in ISI-MIP (20). 
Simulations are binned in temperature bins at AGMT = 1 °C, 2 °C, 3 °C, and 
4 °C (±0.5 °C). For GIM-GCM combinations where the threshold has not been 
crossed by AGMT = 4.5 °C (the highest temperature bin achieved by GCMs in 
this study), a value of 5 is assigned. Consequently, cells with a median sec- 
toral crossing temperature above 4.5 °C are not included in the analysis, 
effectively excluding cells with less than 50% agreement of GIM-GCM 
combinations on the crossing of the respective threshold. See SI Text for 
more details on the sensitivities and uncertainties of the method. If a grid 
cell is identified as having crossed the threshold, the whole area of the grid- 


when the projected values for 2084 are used (Fig. S8). Impacts were simu- 
lated on terrestrial pixels of a global 0.5° mesh (roughly 55 km wide at the 
equator). For an overview of the GIMs used in the analysis, see Tables S1-S4, 
accompanied by a brief discussion of model differences contributing to the 
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cell is assumed to be affected. This process neglects, for example, the sep- 
aration of agricultural and natural vegetation areas in a grid-cell, which is 
below the resolution of the analysis. The spread across GIMs is calculated by 
taking the median over all GCMs for each GIM. The corresponding pro- 
cedure is used for GCMs. 
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